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Abstract 

We calculate total and differential cross sections for J/f; photoproduction in ultrarelativistic 
lead-lead collisions at the LHC energy y/SNN = 2.76 TeV. In the present approach we use a simple 
model based on vector dominance picture and multiple scattering of the hadronic (cc) state in a cold 
nucleus as an example. In our analysis we use both the classical mechanics and quantum (Glauber) 
formulae for calculating atot{J/'fPb) which is a building block of our model. We compare our UPC 
results with ALICE and CMS data. For semi-central collisions (b < Ra + Rb) a modification of 
the photon flux seems necessary. We discuss different physics motivated approximations. We try 
to estimate the cross sections for different centrality bins and for J/tp mesons emitted in forward 
rapidity range (2.5 < y < 4) corresponding to recent ALICE experimental results. Reasonable 
results are obtained but open questions are discussed. 
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I. INTRODUCTION 


Last years the ALICE Collaboration studied production oi J/ip mesons mostly in cen¬ 
tral 208 pp_ 208 pp collisions (see e.g. 0 ) at y/sNN = 2.76 TeV i.e. at the highest avail¬ 
able nucleus-nucleus center-of-mass enereies. At these high energies different theoretical 
mechanisms come into game (see e.g. [2|-l6| and references therein). For instance in Ref. 
jsl the authors emphasized thermal aspect of charmonia production. In order to describe 
detailed distributions in transverse momentum, rapidity or multiplicity several dynamical 
mechanisms have to be taken into account. In contrast to early expectations there are 
mechanisms which lead to both suppression and regeneration oi Jjip quarkonia in cold and 
hot medium. Experimental studies discuss usually so-called nuclear modification factor as a 
function of multiplicity, transverse momentum and rapidity oi J/ip and smooth dependences 
have been observed (a nice summary was presented e.g. in [^). 

In last years there was also an interest in calculating cross sections for exclusive production 
of J/ip mesons in ultrarelativistic heavy ion collisions [9-14|. The main reason for these 
investigations was a better understanding of interaction of a small size (cc) perturbative 
dipole with a nucleus or more explicitly with the nuclear gluons incorporated in different 
ways in different models. In general, this may be related to gluon saturation phenomenon 
which in the case of nucleus may be more efficient than for the case of nucleon due to a 
coherent action of many nucleons (larger gluon occupation). All the previous investigations 
were related to purely ultraperipheral collisions (UPC), i.e. when the nuclei stay intact 
In theoretical calculations UPC means in practice b > Ra + Rb (sum of nuclear radii). All 
the calculations must be therefore performed in the impact parameter space to include this 
condition. 

Can the photoproduction process (photon emission, its fluctuation into cc and subsequent 
rescattering of the dipole (or hadron) in the nucleus) be active also for less peripheral 
collisions i.e. when nuclei collide and break apart producing, at these high energies, quark- 
gluon plasma? A recent analysis of the ALICE Collaboration [l^ presents differential studies 
of inclusive production of J/ip in Pb-Pb collisions exclusively at forward rapidities. At 
low transverse momenta and small multiplicity they observe an enhancement of Raa > 1- 
According to our knowledge this enhancement was not explained so far in the literature. 
We wish to address this issue in the following paper. We will try to discuss the issue in the 
equivalent photon approximation (EPA). In the present first estimate of the cross section 
we will discuss several approximations how to calculate relevant photon fluxes to include 
approximately the physical conditions relevant in the photoproduction of J/ip. This is very 
important for the peripheral and semi-central collisions. For simplicity in the following we 
will use a simple vector dominance model (VDM) combined with multiple scattering of the 
hadronic photon fluctuation for calculating the yA cross section which, as will be shown in 
this paper, describes relatively well the cross section for production oi J/ip mesons in UPCs. 

A very recent analysis of the ALICE Collaboration jl^ conhrms the presence of the 
enhancement at very small J/ip transverse momentum and tries to extract the new not fully 
understood contribution for different centrality bins. We shall try to describe the ALICE 
preliminary data 


16| assuming the mentioned above photoproduction mechanism. 


^ An emission of extra neutrons caused by additional purely electromagnetic interactions could be included 
too. 
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II. DESCRIPTION OF OUR SIMPLE MODEL 


Nuclear photoproduction of a single vector meson V can be understood as a photon fluctu¬ 
ation into hadronic (virtual V meson or quark-antiquark pair) component and its subsequent 
propagation through the second (cold) nucleus and a transformation (fragmentation) into 
the on-shell V meson. 



FIG. 1: Schematic diagrams for the single vector meson production by photoproduction photon- 
Pomeron (left) or Pomeron-photon (right) fusion. Here the word Pomeron and its symbol is an 
abbreviation for multiple diffractive scattering of the hadronic system in the nuclear medium. 


Fig. [T] illustrates a mechanism of a single vector meson production in ultraperipheral 
ultrarelativistic heavy-ion collisions. The cross section for this mechanism is usually written 
differentially in the impact parameter b and in the vector meson rapidity y 

daA iA 2 ^AiA 2 V _ {b, y) dPjp^ {b, y) 

d^dy dy dy 

Above P^^iy^b) or P^p^iy^b) is the probability density for producing a vector meson V 
at rapidity y for hxed impact parameter b of the heavy-ion collision. Probability density 
expresses two-different possibilities of the production of vector meson shown in Fig. [1] Each 
of the probabilities is the convolution of the cross section for 71 A 2 —)■ J/'il)A 2 or 72 Ai 
J/tljAi (the photon emitted from hrst or second nucleus {ujiy = mv '/2 exp(±|/))) and a 
corresponding flux of equivalent photons: 

= , ( 2 . 2 ) 

where N {oJiy-, b) is usually a function of impact parameter between heavy ions ( 6 ) and 
not of photon-nucleus impact parameter. Finally S{b) is an impact parameter dependent 
survival factor which with good precision can be approximated as S{b) ^ 6{b — Ra — Rb)- 
Now we wish to proceed to semi-central collisions i.e. to the case of 6 < Ra + Rb- Then 
the survival factor has to be omitted as we are interested also in situations when colliding 
nuclei break apart. In the present approach we will modify only photon fluxes and leave the 
coherent yA —)• J/'ipA cross section unmodihed. This seems a crude approximation which 
should be reasonable only for rather peripheral collisions b ~ Ra + Rb- We shall try to see 
how far this approximation can be extrapolated down to smaller b. The effective photon flux 
which includes the geometrical aspects can be formally expressed through the real photon 
flux of one of the nuclei and effective strength for the interaction of the photon with the 
second nucleus 

Nd'^{uji,b) = y'iV(a;i,6i) , (2.3) 
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where bi = b + b2. The extra 9{Ra — ^ 2 ) factor ensures collision when the photon hits the 
nucleus-medium. For the photon flux in the second nucleus one needs to replace 1—)-2 (and 
2—)-l). For large b ^ Ra + Rb- (wi, b) ~ iV(a;i, b). For small impact parameters this ap¬ 
proximation is, however, not sufficient. This has some consequences also for ultraperipheral 
collisions, which will be discussed somewhat later in this section. 

Since it is not completely clear what happens in the region of overlapping nuclear densities 
we suggest another approximation which may be considered rather as lower limit. In this 
approximation we integrate the photon flux of the first (emitter) nucleus only over this 
part of the second (medium) nucleus which does not collide with the nucleus-emitter (some 
extra absorption may be expected in the tube of overlapping nuclei). This may decrease the 
cross section for more central collisions. In particular, for the impact parameter 6 = 0 the 
resulting vector meson production cross section will fully disappear by the construction. In 
the above approximation the photon flux can be written as: 


7V(2)(a;,,6)= / iV(a;i, 6 i) 


eiRA - 62 ) X 0(6i - Ra) 


7tR\ 


d%i 


(2.4) 


We shall use the three different approximations how to calculate photon fluxes and three 
different form factors which are main ingredients of the photon flux. We shall start with 


Fpi{g) = 1 


(2.5) 


corresponding to point-like charge. In this approximation the flux of the photons is given 
by the simple formula 

deader Tr^car^ ^ 

where Z is the nuclear charge, u is energy of the photon, r is a distance in the impact 
parameter space between the photon and the emitting nucleus, Ki is a modified Bessel 
function and X = ru/'y. The photon flux integrated over all r is approximately equal to 
the photon flux in the region of r > Ra + Rb- This can be calculated analytically: 



diV(a;) 

do; 


(xKo(x)K,(x) - y [A'?(x) - 


(2.7) 


where y = 2Rauj 

In our calculations we use also a form factor which is Fourier transform of the charge 
distribution in nucleus. Two-parameter Fermi distribution (called equivalently Woods-Saxon 
distribution) 


p{r) 


Po 

l + exp(^) ’ 


( 2 . 8 ) 


where the normalization constant po for the lead nuclei equals to ^0.1572 fm“^ and c = 6.62 
fm, a = 0.546 fm is used in the calculation. The form factor (called here realistic form factor 
for brevity) is calculated then as 


Freal{q^) 


An 

q 


p{r) sin(gr)rdr . 


In the literature often a monopole form factor is used 


F (a^ 

^ mon\H 


X 
X + 


(2.9) 


( 2 . 10 ) 
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which simplifies the calculations, 
square electric radius 


The value of A can be expressed through the root mean 


\/< > 



( 2 . 11 ) 


giving A = 88 MeV for \/< > = 5.5 fm for In our calculation we use the following 

generic formula for calculating the photon flux (see e.g. Ref. [l^) for any nuclear form 
factor 



( 2 . 12 ) 


We shall compare the cross sections for the J/ijj photoproduction using the point-like, 
monopole and realistic form factors (see Table [I]). But let us have first a closer look at the 
photon fluxes. The general situation in the impact parameter space is shown in Fig. |2l 
The two-dimensional vectors b, bi and &2 are distances between colliding nuclei, between 
the photon position and the middle of the first (emitter) nucleus and between the photon 
position and the middle of the second (medium) nucleus, respectively. The production of the 
J/'ip mesons may occur provided the photon (hadronic fluctuation) hits the second nucleus, 
otherwise the (photo)production is not possible. The hatched area of overlapping nuclei 
in the right panel (semi-central collision) of the figure represents the area in the impact 
parameter space for which the situation is not so clear. As will be discussed below we shall 
exclude this region to get lower limit for the cross section for the A A J/'ip reaction. In 
this region of the impact parameter space the quark-gluon plasma is created and its role in 
damping the J/il) production is not clear. The lower limit is obtained by assuming the full 
damping. 



FIG. 2: Impact parameter picture of the collision and the production of the J/'0 meson for ultra¬ 
peripheral (left panel) and for semi-central (right panel) collisions. It is assumed here that the first 
nucleus is the emitter of the photon which rescatters then in the second nucleus being a rescattering 
medium. 

Let us start our discussion by showing standard fluxes used routinely in ultraperipheral 
collisions. Fig. E] presents the two-dimensional photon fluxes (see Eq. fl2.12p l as a function 
of the distance between two colliding nuclei (6) and the energy of the photon (cj). One 
can observe that the difference between photon fluxes obtained with realistic (left panel) 
(or monopole (middle panel)) form factors and the result for point-like photon source (right 
panel) is huge for b < Ra (especially for 5 0). 

Fig. m presents the realistic photon flux which is calculated including extra absorption 
effects in three different ways. The left panel shows the result which is obtained from 
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FIG. 3: Standard photon fluxes calculated for realistic (left panel) monopole (middle panel) and 
point-like (right panel) form factors. 



FIG. 4: Two-dimensional distributions of the photon flux in the impact parameter b and in the 
energy of photon u for three different conditions (more in the text). 


Eq. (12.3p . This limit allows a production of vector mesons only inside the second nucleus 
(the one which does not emit the photon). The second panel is obtained with a condition 
which follows from the dehnition of the absorption factor. 6{bi — Ra) in Eq. (12.4p allows a 
production of the meson everywhere except of the inside of the tube dehned by the emitter. 
The last panel shows the result which is obtained from Eq. (12.4p . Here we have included 
both conditions applied separately in the two previous panels of the hgure. This condition 
can be understood as follows. For centrality smaller than 100% we consider only these cases 
when particle can be produced inside the medium excluding this part of the medium which 
coincides with the emitter in the impact parameter space. 

We think that it is more pedagogical to inspect the ratio of the effective flux in the new 
approach to the standard photon flux (called below N^^'>{u,b)). Fig. E] presents 

the ratio of the ’’new” to the standard flux for three types of the nuclear form factor. More 
precisely, it is a ratio of the result presented in the last panel of Fig. |l]to the result shown 
in the hrst panel of Fig. |3l The three panels represent cases with realistic (left panel), 
monopole (middle panel) and point-like (right panel) form factors. Here we consider lead- 
lead collisions at the LHC energy ^/snn = 2.76 TeV (7 = 1471). It is shown that the ratio 
of b) to b) is almost 1 in the region of large b for all form factors. However, it 

is not the case for the semi-central collisions. These differences very weakly depend on the 
energy of the emitted photon. Thus interesting seems to be the cross section for different 
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FIG. 5: The ratio of the differential photon fluxes in the impact parameter b and energy of the 
photon uj. The left panel shows the case with the realistic form factor, the middle panel monopole 
form factor and the right panel the form factor for point-like charge. 


ranges of centrality. 

Centrality is a parameter which is often nsed to characterize the collision of nnclei at 
high energies. To efficiently compare onr resnlts with the preliminary ALICE data for semi¬ 
central collisions [l^ we need a simple relation between impact parameter and the centrality. 
In Ref. [l^ the anthors gave an approximate and practical geometric relation between them. 
In general, centrality of collisions depends on the impact parameter b and a total inelastic 
nnclens-nnclens cross section 

7r62 (AT) 

^inel 

In this case the centrality corresponds to a mnltiplicity higher than N. b{N) expresses the 
valne of the impact parameter for which the average mnltiplicity fnlhlls the dependence 
n{b) = N. In a pnrely geometrical pictnre {(Tinei = vr we get: 


(2.13) 



c = 


62 

■ 


(2.14) 



FIG. 6: Gentrality of nuclear collisions as a function of the impact parameter as obtained from 
Eq. (I2TTD . 
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Fig. |6] presents the centrality as a fnnction of impact parameter. The large value of the 
centrality corresponds to the large value of the impact parameter. The larger value of b the 
more peripheral collisions. 

Let us concentrate now for a while on the second ingredient of the model (see Eq. fl2.2p l 
- the (J^A^j/ipA cross section. In our present calculations we use the following sequence of 
equations: 


d(T {'jp —)■ J/tjjp] t = 0) 
dt 


— bj/^ijXj/^W^p 




da{J/'il)p^ J/ipp]t = d) fj/ijj da (yp-)■ J/V’p; t = 0) 


dt 


= IStt 


Ta (r) = 


dvrOem dt 

da {J/'ipp J/i’P'i t = 0) 
dt 

[ p ^\/r2 + dz , 


da (yA —)■ J/iIjA] t = 0) 
dt 




4/ 


(2.15) 

(2.16) 

(2.17) 

(2.18) 

(2.19) 


JIP 


Constants for the yp —>■ J/'ijjp production are obtained from a £t to HERA data (Ref. 
fl^): bj/^ = 4 GeV“^, Xj/^ = 0.0015 ph, ej/p = 0.8 and vector-meson coupling square is 
= dvr • 10.4. 

The calculation of atot{J/'ip^) requires incorporating multiple scattering of J/'ip in the 
nuclear medium. In the present paper, the calculations are done using both classical 
mechanics {CM) and quantum mechanical {QM) Glauber formula for atot{J/4’A) cross 
section. ”A quantitative” and spectacular difference between these approaches can be 
shortly summarized as follows. In the black disk limit the classical mechanics approach 
leads to the total cross section equal to and the quantum mechanical approach im¬ 


plies atot{J/'4’A) = 27rR\ j20|. The classical and quantum mechanics expressions for the 


O'tot (J/i’A) cross section read: 


i'J/i^A) = j d^r (1 - exp {-atot (J/i^p) Ta (r))) , (2.20) 

{J/i’A) = 2 j d^r(^l- exp {Jl^p) Ta (r)^ ^ , (2.21) 

respectively, where r is the distance in the impact parameter space of the photon (or cc 
fluctuation) from the middle of the nucleus-medium. These formulae are used to calculate 
the yA —)■ J/ijjA cross section, the main ingredient of our whole approach. 

Finally, the total cross section for yA —)■ Jji\}A reaction can be written as 


fljjA 


dO'-yA^JllpA {t 0 ) 

dt 


tmax 

J dt\FA{t)f . 


— OO 


( 2 . 22 ) 


This cross section is actually a function of energy in the yA system. The factor \FA{t)\ 
appears here due to assumed coherent rescattering of the qq dipole off a nucleus. A very 
good approximation is to use the realistic nuclear charge form factor which is dehned in 
Eq. fl2.9p . The squared four-momentum transfer: t = —q^ = —{^jfpl {‘^^lab))'^- 














III. NUCLEAR CROSS SECTIONS FOR UPC AND MORE CENTRAL COLLI¬ 
SIONS 


We start verification of our model for purely ultraperipheral collisions. Fig. [7] presents 
results for coherent photoproduction in the 208 pp_|_ 208 p]^ UPC at = 2.76 TeV. We 

show our results for classical and quantal rescattering in the nucleus-medium and for different 
nuclear form factors. The result with the monopole form factor strongly overestimates 
the ALICE and CMS data. The result with the quantal rescattering is about 15% larger 
than that for the classical rescattering. The difference is much smaller here than for the 


photoproduction of meson [2l|]. The CMS ’’data point” was obtained by correcting a 


real data point with at least one neutron in zero-degree calorimeter (ZDC) using the Monte 


Carlo program STARLIGHT [22(]. The ALICE data points are taken from Refs. j2J, [25 


and the CMS data point is from Ref. [23|. Relatively good agreement is obtained for the 


forward rapidity region provided the realistic nuclear form factor is used which supports 
application of the model also for more central {b < Ra + Rb) collisions especially in the 
considered rapidity range. 



FIG. 7: Differential cross section for coherent production of J/'ijj meson in UPC as a function of 
rapidity of the J/ij; meson compared with the ALICE and CMS data points. We show results for 
both realistic and monopole form factor, each of them used consistently in Eq. (j2.12p and (I2.22|) . 

In Fig. |8]we show the nuclear cross section as a function of the impact parameter also 
for its small values {h < Ra + Rb) he. for the semi-central collisions. The different lines 
correspond to different approximations of photon fluxes within our approach as described 
in the figure caption. The dashed and solid lines represent upper and lower limit for the 
cross section. At larger values of impact parameter h the cross sections obtained with the 
different fluxes practically coincide. At 6 < Ra + Rb fhe different approximations give quite 
different results. The standard approach used in the literature for UPC (Eq. fl2.12p l when 
naively applied to the semi-central collisions overestimates the cross section. We will return 
to this issue in the following. 

In Tab. m we have collected the cross sections for six ranges of collision centrality used in 
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FIG. 8 : Differential cross section for photoproduction of J/ij: meson as a function of impact pa¬ 
rameter for y/sNN = 2.76 TeV. Different lines correspond to different approximations: dotted 
- standard UPC approach, dashed - first approximation/correction (upper limit), solid - second 
approximation/correction (lower limit). Here realistic (charge) nuclear form factor was used. 


Centrality range [%] 

0-10 

10-30 

30-50 

50-70 

70-90 

0-100 

UPC 

^REAL 

1.19 

3.36 

2.58 

1.57 

0.97 

10.08 

16.32 

[mb] 

3.64 

3.99 

2.16 

1.36 

0.99 

12.58 

16.26 

[mb] 

211.00 

6.11 

2.46 

1.50 

1.06 

222.68 

16.43 

dcr/dyj/^ [pb] {Ay = 1.5) 

ALICE data 

<318 

<290 

73 

±44^26 

58 

±i 6 ;«o 

59 

4-11+7 



[pb] {Ay = 1.5) 

111 

282 

196 

105 

55 



do-M^^^^^/dyj/^ [pb] {Ay = 1.5) 

398 

358 

158 

86 

56 



da^^/dyj/^ [pb] {Ay = 1.5) 

9984 

579 

185 

92 

58 




TABLE I: Integrated and differential cross section for the photoproduction of J/'ijj in Pb + Pb 
collisions for = 2.76 TeV calculated with a simple flux factor according to Eq. (I2.12p for the 

case when realistic (REAL), monopole (MONOP) and point-like (PL) form factors are used. We 
show results for different ranges of centrality as well as for ultraperipheral events. The differential 
cross section is compared with the ALICE experimental data [1^. Ay = 1.5 means y G (2.5,4.0) 
for brevity. 


16l | as well as for UPC. We show total (full range of rapidity) and differential (da/dyj/^) 


cross section for rapidity range of the ALICE experimental studies (i/j/p = (2.5—4)). Here we 
take the photon flux given by Eq. fl2.12p . For each case, the largest cross section appears for 
centrality range 10 — 30% which corresponds to 6 ~ (4.4—7.6) fm according to formula fl2.14p . 
One should be aware that the results for the point-like form factor {F{q) = 1) at very small 
centrality are not physical. Here there is a very strong dependence on impact parameter 
b. The photon flux for very small value of b tends to inhnity. Looking at the agreement 
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of our results with the experimental ALICE data, one can observe that calculations with 
realistic form factor are in fairly good agreement with the experimental data. For centrality 
larger than 30% the cross section obtained with the monopole form factor is closer to the 
experimental value. However, for each form factor the cross section is larger than that given 
by the ALICE Collaboration. Only for centrality range 70 — 90% we get a good agreement 
with the experimental result. 

It is worth to note here that the nuclear form factor actually appears in two formulas 
in the calculation of the cross section for the AA AAJ/ip reaction. A Erst one is for 
the flux of photons (see Eq. fl2.12p i and the second one is for the cross section for the 
7 % —)■ J/ijjA reaction (see Eq. fl2.22p i which is integrated over squared four-momentum 
transfer. Only hard sphere form factor was used so far in the literature in the context of 
the J/ip photoproduction [^. The results presented in Table |T] are calculated for the case 
when realistic charge distribution is included both in Eq. fl2.12p and in fl2.22p . and for the 
case when the monopole or point-like form factors are used only in calculating the photon 
flux and realistic form factor is used in Eq. fl 2 . 22 p . We suppose that a consistent calculation 
should use the same form factor in both places (Eqs. 02 . 121 ) (or 02.41) / 02.3p l and 02 . 22 p ). 


Centrality range [%] 

0-10 

10-30 

30-50 

50-70 

70-90 

0-100 

UPC 

|mb| 

8.23 

9.03 

4.89 

3.08 

2.45 

28.50 

38.42 

[mb] 

17 037.50 

926.26 

379.30 

223.40 

158.50 

18795.50 

2838.30 

du/dyj/^ [/ib] (Ay = 1.5) 

ALICE data 

<318 

<290 

73 

±44^26 

58 

±16l?o 

59 

4-11+'^ 



^^MUNUF [/ib] (Ay = 1.5) 

892 

817 

365 

201 

132 



dcr^'^/dyj/^ [/ib] (Ay = 1.5) 

1 488 700 

88 230 

28 800 

14 660 

9 370 




TABLE II: Integrated and differential cross section for the production of J/'ip in Pb + Pb collisions 
for y/sjsTN = 2.76 TeV for the case when monopole (MONOP) and point-like (PL) form factor is 
included in photon fluxes and in the total cross section for the yA —>• J/'ij^A reaction. We show 
results for different ranges of centrality as well as for ultraperipheral case. The differential cross 
section is compared with the ALICE experimental data 
brevity. 


16l |. Ay = 1.5 means y G (2.5,4.0) for 


In Table im we have collected the nuclear cross section for the PbPb —)■ PbPbJ/ip reaction. 
This table is very similar to the previous one, but now we use the monopole and point-like 
form factor also in the formula which describes the total cross section for the 7 A —)■ J/xpA 
reaction fl2.22p . These numbers are much larger than in the previous approach. We conclude 
that the point-like approximation cannot be used to calculate the 7 P 6 —)■ J/ipPb cross 
section because it leads to unphysical behaviour for b < Ra + Rb- 

Table IIIII and Table nYi collect the cross sections for different centrality bins for the hrst 
(Eq. fl2.3p ) and second (Eq. fl2.4|) ) approximation of the photon flux, respectively. A fairly 
good agreement with the ALICE data is obtained with the second approximation for all 
measured centrality bins. 

We summarize the numbers from the tables in Fig. [9l We present both statistical and 
systematic error bars (shaded area). The ALICE Collaboration could not extract actual 
values of the cross section for the two lowest centrality bins. The results for standard 
photon flux exceed the ALICE data. Rather good agreement with the data is achieved for 
the photon flux obtained with the realistic nucleus form factor. The fact that this 
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Centrality range [%] 

0-10 

10-30 

30-50 

50-70 

70-90 

0-100 

^REAL 

1.42 

2.47 

2.08 

1.60 

1.24 

9.37 

dn/dyj/^ [pb] {Ay = 1.5) 
ALICE data 

<318 

<290 

73 

±44^2? 

58 

±16l?o 

59 

4-11+7 


do-^^'^^/dyj/^ [pb] {Ay = 1.5) 

123 

201 

160 

116 

84 



TABLE III: Integrated and differential cross section for the production of in Pb+Pb collisions 
for y^SNN = 2.76 TeV calculated with the help of the first approximation of the photon flux 
N^^\u},b) (Eq. (I2.3p i for the realistic form factor. The differential cross section is compared with 
the ALICE experimental data 1^. Ay = 1.5 means y G (2. 5,4.0) for brevity. 


Centrality range [%] 

0-10 

10-30 

30-50 

50-70 

70-90 

0-100 


0.38 

1.19 

1.29 

1.21 

1.11 

5.47 

dn/dyj/^ [pb] {Ay = 1.5) 
ALICE data 

<318 

<290 

73 

±441P 

58 

±16l?o 

59 

4-1 1 +7 


^^REAL [pb] {Ay = 1.5) 

30 

88 

91 

82 

72 



TABLE IV: Integrated and differential cross section for the production of J/ij) in Pb-\-Pb collisions 
for y/SMN = 2.76 TeV calculated with the help of the second approximation of the photon flux 
N^'^\u},b) (Eq. (I2.4p i for the realistic form factor. The differential cross section is compared with 
the ALICE experimental data 1^. Ay = 1.5 means y G (2. 5, 4.0) for brevity. 



EIG. 9: da/dy cross sections for different centrality bins. Theoretical results for different models of 
the photon flux are compared with the preliminary ALICE data 1^. The shaded area represents 
the experimental uncertainties. 


lower limit exceeds somewhat the ALICE data may be due to the fact that the coherent 
7 A cross section was used in our EPA calculation whereas only spectators may be 

active in the production of J/H/. 
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IV. CONCLUSIONS 


So far photoproduction of vector mesons was considered only in the context of ultrape¬ 
ripheral collisions, i.e. in the case when nuclei survive the collision and appear in the hnal 
state, although the exclusivity is not checked in practice and experimentally a limited rapid¬ 
ity gap condition is imposed only. In this paper we have presented a hrst theoretical study 
of photoproduction mechanism also in the case when nuclei collide and produce quark-gluon 
plasma and as a consequence considerable number of hadrons, mostly pions is produced. On 
theoretical side, the nuclear photoproduction in UPC is treated in the equivalent photon ap¬ 
proximation with photon fluxes and photon-nucleus cross section being the basic ingredients 
of the approach. 

How to calculate contribution of the photoproduction mechanism for peripheral and semi¬ 
central collisions (6 < Ra + Rb) is not completely clear. We assume that the whole nucleus 
produces photons. The photon (or hadronic photon fluctuation) must hit the other nucleus 
to produce the J/ijj meson. The question arises how to treat the region of overlapping 
colliding nuclei in the impact parameter space where some absorption may be expected. We 
decided to include the effect of the ” absorption” by modifying effective photon fluxes in the 
impact parameter space by imposing several additional geometrical conditions on impact 
parameters (between photon and nuclei and/or between colliding nuclei). 

In the present paper, as an example, we have considered a vector-dominance based model 
which includes multiple scattering effects in addition. Any other model/approach can be 
applied in the future. 

First we have compared results of our calculations with existing ALICE and CMS exper¬ 
imental data for UPC. We have obtained a reasonable agreement with the data especially in 
the forward rapidity region when realistic nuclear form factor is used. We have calculated 
the dependence of the total and differential cross section on the impact parameter (the dis¬ 
tance between colliding nuclei in the plane perpendicular to the collision axis). By correcting 
standard photon fluxes valid only for UPC by collision geometry we have calculated cross 
section for different centrality bins adequately for the recent ALICE Collaboration analysis. 
Our results have been compared with their preliminary data. We have obtained a reason¬ 
able agreement for peripheral and semi-central collisions and set a lower and upper limits 
for the cross section for the semi-central collisions within our approach. Our lower limit is, 
however, model dependent. Since in our calculations we have used coherent 7 A —)■ Jji\}A 
cross section our lower limit may be overestimated especially for small impact parameters 
where only so-called spectators are active for the J/'ip production. The time picture of the 
whole process is not clear to us in the moment. The rather reasonable agreement of our 
quite simplihed approach with the preliminary ALICE data could suggest that the ’’coher¬ 
ent” (assumed by the formula used for the 7 A —)■ J j ip A process) scattering of the hadronic 
fluctuation happens before the nucleus undergoes the process of deterioration due to nucleus- 
nucleus collision and before the quark-gluon plasma is created. If this is not the case the 
’’coherent” scattering must be replaced by incoherent one. Then the transition from UPC 
to semi-central collisions as a function of impact parameter would be more abrupt then in 
the present approach. 

In the present paper we have performed analysis for a forward rapidity range. There 
the J/ip quarkonia are emitted forward with large velocity therefore they could potentially 
escape from being melted by quark-gluon plasma. At midrapidities the situation can in 
principle be different. It would be therefore valuable if the ALICE Collaboration would 
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repeat their analysis also in the midrapidity range and verify the pt ^ 0 enhancement. 

The photoproduction of J/V’ seems rather special compared to photoproduction of other 

(jg-c/pc (y=0) 

objects. Here the corresponding cross section (— pbPb^pbj^bj/^ - ^ 2 j^ib, see Fig. [7]) is 


of similar order of magnitude as the one for J/ijj produced in individual nucleon-nucleon 
collisions or produced in quark-gluon plasma (simple estimate including ~ 4 


dy 


/ib (Ref. [3]), estimated number of binary collisions and experimental nuclear modihcation 
factor (Ref. [I|) gives ^ ^ fQj, minimum-bias collisions). In our simple 


^^photoprod. _ Q ^ 

estimate we get —^ (0.91 — 1.42) mb, where the two numbers are 

obtained with and fluxes, respectively. These numbers are of similar size as the 
contribution of other mechanisms discussed above. For (photo)production of other objects 
the situation may be different. This was not discussed in the literature so far and requires 
dedicated studies. 
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